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^ ■ ABSTRACT 

rr . We report 2.1 mm observations of the Sunyaev-Zel'dovich (S-Z) effect; these obser- 

^_. I vations confirm our previous detection of a decrement in the Cosmic Microwave Back- 

^ ' ground intensity towards the cluster Abell 2163. The S-Z data are analyzed using the 



relativistically correct expression for the Comptonization. We begin by assuming the 
intracluster (IC) gas to be isothermal at the emission weighted average temperature de- 
tJ" I termined by a combined analysis of the ASCA and GINGA X-ray satellite observations. 

^ ■ The results of ROSAT/PSPC observations are used to determine an isothermal model 

^ ■ for the S-Z surface brightness. Fitting to this model, we determine the peak Comp- 

O ! tonization to be yo = 3.73l;gi x 10""^. The uncertainty includes contributions due to 

statistical uncertainty in the detection, instrumental baseline, calibration, and density 
model. Combining the X-ray and S-Z measurements, we determine the Hubble constant 






p i' to be Ho{qo = |) = GOlllJkms ^Mpc ^, where the uncertainty is dominated by the sys- 



O 



tematic difference in the ASCA and GINGA determined IC gas temperatures. ASCA 



?-j I observations suggest the presence of a significant thermal gradient in the IC gas. We 

^ ■ determine Hq as a function of the assumed IC gas thermal structure. Using the ASCA 

determined thermal structure and keeping the emission weighted average temperature 
the same as in the isothermal case, we find i/o(Q'o = 1) ~ 78lilQ]iias~^Mpc~^ . Including 
rS I additional uncertainties due to cluster asphericity, peculiar velocity, IC gas clumping, 

c3 ■ and astrophysical confusion, we find i^ol'Zo = 2) ~ 78^40 kms~^Mpc^^. 
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Introduction 



the Einstein, GINGA, ROSAT, and ASCA satellites 



Conipton scattering of rhe cosmic microwave back- 



ground (CMB) radiation by hot inrraclusrer (IC) gas 
- the Sunyaev-Zel'dovich (S-Z) effect - gives rise to 
an observable distortion of the CMB spectrum which 
can be used as a sensitive cosmological probe. The in- 
tensity change caused by the scattering has a thermal 
component, due to the random motions of the scat- 
tering electrons ( Sunyaev fc Zcl'dovich 1972 ), and a 
kinematic component, due to the bulk peculiar ve- 



locity of the IC gas ( Sunyaev fc Zel'dovich 



When combined with X-ray measurements, the am- 
plitude of the thermal component can be used to de- 



termine the Hubble constant (Gunn 1978; Cavaliere 



et al. 1979). The ratio of the amplitudes of thermal 



and kinematic components can be combined with the 
X-ray temperature to determine the radial component 
of the cluster peculiar velocity (^unyaev fc Zel'dovich 



1980; Rephaeli & Lahav 1991). Because the surface 



brightness of the S-Z effect is independent of redshift, 
this method for determining the Hubble constant and 
cluster peculiar velocities can be applied to distant 
clusters, provided the requisite X-ray data can be ob- 
tained. 

The thermal and kinematic components of the S-Z 
effect have distinct spectral shapes and can be sepa- 
rated from one another by measurements at two or 
more millimeter (mm) wavelengths. Astrophysical 
confusion for the S-Z effect exhibi ts a broad mini- 
mum at wavelengths of 2 — 3 mm ( Fischer fc Lange 
1993 ). Several significant detections of the S-Z effect 



at centimeter (cm) wavelengths have been reported 
in recent years. Some of these measurements have 
been combined with X-ray observations to compute 
values for Hq (Birkinshaw, Hughes fc Arnaud 1991 



Birkinshaw fc Hughes 1994; Herbig et al. 1995; Jones 



et al. 1993; Myers et al. 1995; for a recent review 



see Rephaeli 1995a). The prospect of separating the 
thermal and kinematic components of the S-Z effect 
while minimizing astrophysical confusion makes ob- 
servations at mm wavelengths attractive. 

We have developed a novel instrument and observ- 
ing strategy with high sensitivity to the S-Z effect in 
several millimeter bands (Holzapfcl et al. 1996a). Us- 
ing this instrument, we have detected the S-Z effect 
in the direction of the cluster Abell 2163 (Wilbanks 
et al. 1994| - hereafter W94). A2163 is a ri ch, moder- 
ately distant (z = .201, poucail et al. 1995 ) cluster of 
galaxies, which has been observed in the X-ray with 



(Arnaud et al. 1992; Elbaz et al. 1995; Markevitch 



et al. 1996). Its X-ray luminosity and temperature 



are among the highest of any cluster yet studied, im- 
plying a high IC gas density and temperature, and 
correspondingly high S-Z surface brightness. The low 
declination and modest angular extent of this clus- 
ter make it a good candidate for observations with 
our instrument, which operates by drift scanning (see 
Sections f4.1.| fc [4.2.[ ). In this paper, we combine ob- 
servations of the S-Z effect and X-ray emission to de- 
termine Hq. This is the first determination of Hq 
from millimeter wavelength measurements of the S-Z 
effect. 

The determination of Hq presented here includes 
several refinements with respect to previous work. 

1.) We use a relativistically correct calculation of 
CMB Comptonization in clusters. The nonrelativis- 
tic treatment is inaccurate at high frequencies and 
high IC gas temperatures (see Section 2T.). Use of 



the relativistically correct formulae in the analysis of 
2.1 mm measurements of A2163, the hottest known X- 
ray cluster, is essential in order to derive an accurate 
value of Hq . 

2.) We use the combination of the ROSAT derived 
density profile and the S-Z surface brightness to con- 
strain models for the large scale thermal structure of 
the IC gas. The allowed range of thermal structure 
models is compared to that determined from the com- 
bination of ASCA and ROSAT measurements. We 
have determined j/o a-nd Hq as a function of the cluster 
thermal structure, specifically for the thermal struc- 
ture measured by ASCA. 

3.) Finally, we include a detailed analysis of addi- 
tional contributions to the uncertainty in Hq due to 
uncertainties in the distribution and thermal struc- 
ture of the IC gas, the possibility that the cluster 
has an appreciable peculiar velocity, and astrophysi- 
cal confusion to the S-Z effect. 

In Section E], we briefly review the relativistic 
Comptonization calculation and outline the method 
used to determine Hq from S-Z and X-ray measure- 
ments. The analysis of the X-ray data is outlined 
in Section p\. The SuZIE instrument, its calibration, 
and observations of A2163 are described in n\. In 
Section pi we describe the analysis of the S-Z data 
under the assumption of isothermal gas. We combine 
the isothermal S-Z and X-rayresults to determine the 
Hubble constant in Section |6]. In Section ^ we dis- 



cuss the large scale thermal structure of the IC gas 
and its effects on our results. Additional contribu- 
tions to the uncertainty in Hq are discussed in Sec- 
tion pj our conclusions are summarized in Section p]. 

2. Theory 

2.1. Relativistic Comptonization in Clusters 

The Sunyaev & Zel'dovich (1972) treatment of 
CMB Comptonization by hot IC gas is based on a 
solution to the Kompaneets (1957) equation, a non- 
relativistic diffusion approximation to the full kine- 
matic equation for the change of the photon distribu- 
tion. The expression derived for the intensity change 
- an expression which has been used in virtually all 
works on the S-Z effect in clusters - is valid only at 
low gas temperatures. A more accurate calculation 
of the thermal effect requires use of the relativistic 
form for the electron velocity distribution and the full 
expression for the scattering probability, as has been 
presented by Rephaeli (1995b). Use of the correct rel- 
ativistic calculation in this work is essential because of 
the high observing frequency and the unusually high 
IC gas temperature in A2163. 

In the nonrelativistic limit, the CMB intensity 
change due to Compton scattering by a hot electron 
gas (with no bulk velocity) is (Zel'dovich & Sunyaev 
19681) 



AInr = loygix), (1) 

where x = hv/kTi^,, To is the CMB temperature, and 
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The spectrum is given by the function, 
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which va nishes at xp = 3.83 (;^o — 217 GHz) for Tq = 
2.726 K (iMather et al. 1994| ). The Comptonization 
parameter is given by, 



y 



kT, 



n^aTdl , 



(4) 



where Ue and Te are the electron density and tem- 
perature, ut is the Thomson cross section, and the 
integral is over a line of sight through the cluster. 

As has been shown by Wright (1979), the exact 
scattered photon distribution can be calculated using 



the Chandrasekhar (1950) expression for the scatter- 
ing probability and the relativistically correct form for 
the Maxwellian electron velocity distribution. Inte- 
grations over the scattering angle, electron velocities, 
and incident (Planck) spectrum, yield an expression 
for the intensity change (Rephaeli 1995b), 



A/ = /f 



f^ix,Te)dT, 



(5) 



where dr = rieOTdl. The spectrum is given by 



*(.T, T,) 



(e^ - 1) 



[<i>(.T,Te)-l], (6) 



where <I'(x,Te) is a three dimensional integral over 
electron scattering angles, velocities, and the change 
in frequency of the scattered photons which is fully 
specified in Rephaeli (1995b). This calculation is 
exact to first order in r. In Fig. ^, we compare 
the relativistic and nonrelativistic expressions, at the 
GINGA -I- ROSAT derived isothermal gas tempera- 
ture in A2163, kTe = 14.6 keV (|Elbaz et al. 1995|). 
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Fig. 1. — The spectral forms of the intensity change 
due to the thermal Sunyaev-Zel'dovich effect and the 
non-rclativistic approximation. The solid line shows 
A//t, and the dotted line shows A/„r/T (both in 
units of 2[kTQf/{hcf) for kTe = 14.6 keV. 

At the central frequency of the SuZIE 2.1mm 
band {v ~ 142 GHz), the ratio of the brightnesses 



is AI/AI„r = -91. Finally, in order to facilitate com- 
parison of our results with low frequency observations, 
we combine equations and H to express an effective 
Comptonization in terms of the measured intensity 
change, 

^ A/ I^dr 
y /o J^{x,T,)dT- ^'' 

All results for the Comptonization presented in this 
paper are computed using this expression. 

2.2. Hubble Constant from S-Z and X-Ray 
Measurements 

The determination of Hq from the combination of 
X-ray and S-Z data has been treated by Birkinshaw, 
Hughes, & Arnaud (1991) - hereafter BHA. Their 
treatment is adopted here with some exceptions; we 
express the S-Z decrement in intensity rather than in 
Ray leigh- Jeans temperature and use a relativistically 
correct expression for the Comptonization. 

Hot IC gas is detected primarily through its ther- 
mal bremsstrahlung X-ray emission. The X-ray sur- 
face brightness at energy E can be expressed as the 
line of sight integral (BHA equation 3.1), 



1 



Ox 



47r(l + 2)' 



■ f nlA{E,Te)dl., 



(8) 



where Ue is the electron density, A{E, Te) is the spe- 
cific spectral emissivity, and z is the cluster redshift. 
We express the electron temperature, density, and 
spectral emissivity in terms of the product of a ref- 
erence value and a dimensionless form factor (BHA 
equations 3.3 - 3.5): 



Ae{E, T,) ^ AeofAiE, Te) 



(9) 
(10) 

(11) 



where 9 is the angle from a reference line of sight, 
(p is the azimuthal angle about that line of sight, ^ 
is the angular distance along the line of sight, and 
rieo and T^o are the central gas density and temper- 
ature. We extend this treatment to include a term 
describing the temperature dependence of relativistic 
Comptonization, 



*(a;, re) = *o/*(a;, Te 



(12) 



The expressions for the X-ray and S-Z surface bright- 
nesses, equations || and |[ then become: 
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= NxOx , 
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lo'i'oneocrTd.A / d^fnh , 


(15) 


= Nsz^sz , 


(16) 
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where dA is the angular diameter distance. The struc- 
tural information for the cluster is contained in the 
angles: 

0x= I dull A, (17) 



QSZ = j dUnh ■ 

The normalization factors are then defined as: 

neoAeO^A 



Nx 



47r(l-h2) 



Nsz = h^meQcrrdA ■ 



(18) 



(19) 



(20) 



Eliminating the unknown central electron density be- 
tween these equations, we can express dA in terms of 
measurable quantities. 



dA 



Nh 



Ae 



Nx J 47r(l + z)^[/o*oH' 



(21) 



which is an implicitly expressed generalization of 
BHA's equation (3.12) valid at all frequencies and gas 
temperatures. Combining Nx and Nsz with a mea- 
surement of Teo, which appears in the expressions for 
Aeo and ^Oj we can then solve for ifo('Zo) ex: l/rfA- 

Because we have no information about the distri- 
bution of the IC gas along the line of sight, we are 
required to make some simplifying assumptions about 
the temperature and density form factors. We assume 
that the density form factor can be described by a 



convenient empirical model (Cavaliere et al. 1976) 



Ue) = 1 



ei 



(22) 



in which [3 and the cluster angular core radius {6c) 
are left as free parameters. The thermal structure 
of the IC gas is also assumed to be radially symmet- 
ric. Together these assumptions make possible the 



determination of /„ and fx from the observed X-ray 
surface brightness as a function of X-ray energy. 

Because A{E, Tg) and ^(x, Tg) are non-linear 
functions of the temperature, the shapes of the X- 
ray and S-Z surface brightness profiles depend not 
only on the temperature form factor but also on its 
reference value. In general, the form factors cannot 
be determined from monochromatic imaging data if 



the l cluster gas is not isothermal. However, at ener- 
gies much lower than the electron thermal energy (as 
in the ROSAT/PSPC energy band), the X-ray emis- 
sivity of the IC gas depends weakly on the IC gas 
temperature. Therefore, even for a non-isothermal IC 
gas, a good approximation to the density form factor 
can be determined from the X-ray surface brightness 
profile. 

The S-Z surface brightness is roughly proportional 
to the electron pressure integrated along the line of 
sight through the cluster. It is possible to constrain 
spherically symmetric models for fx by integrating 
the product of /„ (determined consistently from the 
X-ray data for each fx model) and f^{x,Te) along 
the line of sight to determine a model for the S-Z sur- 
face brightness, which is then fit to the 2.1mm data. 
In Section [zj, we use this method to constrain mod- 
els for the IC gas thermal structure. Presently, much 
stronger constraints on the thermal structure are de- 
termined from the spatially resolved X-ray spectra 
obtained by the ASCA satellite. 

3. X-ray Analysis 

A2163 was the target of pointed observations 
by the GINGA and ROSAT satellites. A de- 
tailed combined analysis of these observations has 
been published previously ( Elbaz et al. 1995 ). The 
ROSAT/PSPC and GINGA spectra were simultane- 
ously fit to a redshifted isothermal plasma emission 
model where the hydrogen column density along the 
line of sight (A^h), heavy metal abundance, total emis- 
sion measure, and IC gas temperature were free pa- 
rameters. The best fit emission weighted temperature 
was found to be Tg = 14.6lokeV, at 90% confidence 



dElbaz et al. 19951 ). 

The ROSAT/PSPC data were used to determine 
the spatial dependence of the X-ray emission. In 
Fig. g, we show the X-ray surface brightness as de- 
termined by the PSPC. The peak surface brightness 
was found at 161^15^46% -06° 08' 55" (J2000). The 
radi ial profile of the surface brightness (corrected for 



vignetting) was found by summing annuli about this 
central value, with significant emission detected up 
to 18'. The surface brightness profile was fit with the 
combination of an isothermal j3 model (equation p2h 
convolved with the instrument point spread function 
(PSF) and a constant background. The best fit den 
sity model parameters were found to be /3 = .62 



and6>c = 1.20 ±0.075' at 90% confidence (Elbaz et al. 



1995) 
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Fig. 2. — X-ray image of the A2163 cluster ob- 
served with the ROSAT/PSPC and corrected for 
vignetting. The contour levels are logarithmically 
spaced with the peak brightness corresponding to 
9.4 X 10~^ cts arcmin^^. Superimposed are the SuZIE 
drift scans made in 1994; one row of detectors passes 
over the X-ray center and the other 2.2' to the south. 
Scans were alternately begun 12' (shaded) and 18' 
ahead of the X-ray cluster. 

Using the ASCA satellite it is, in principle, possible 
to measure large-scale variations in the temperature 
of the IC gas. Combining the ASCA/SIS-I-GIS spec- 
tra with the ROSAT/PSPC derived density model, 
Markevitch et al. (1996) determined the tempera- 
ture of the IC gas in three radial bins of — 3, 
3 — 6, and 6 — 13 X-ray core radii to be 12.2^|^'2 keV, 
ll.StlgkeV, and 3.8tJ;gkeV at 90% confidence. 

As noted by Markevitch et al. (1996), the 
ASCA determined temperature is systematically 
lower than that found by GINGA. The emission- 
weighted ASCA/GIS-I-SIS average temperature for 
the combination of the 3 radial bins is Tg = 11.2 ± 
1.1 keV at 90% confidence (Markevitch, private com- 
munication), while an analysis of the GINGA data 



gives Te = 14.61 gkeV, also at 90% confidence (Elbaz 



et al. 1995). These two temperature determinations 



are inconsistent at 99.99% confidence. It is beyond 
the scope of this paper to try to understand the ori- 
gin of this discrepancy, part of which may originate 
in the cahbration of the instrument effective areas. 

In Section ^.2j , we estimate the uncertainty in the 
central and emission-weighted average temperatures 
of the IC gas taking this systematic difference into 
account. In addition, we discuss the effect of the un- 
certainty in temperature on the other X-ray parame- 
ters relevant for the determination of the Hubble con- 
stant. In Section 3.1., we consider the ASCA deter- 



mined thermal structure and introduce a simple two 
parameter model for the IC gas temperature profile. 



In Section 3.3.3., we repeat the analysis of the ROSAT 



imaging data consistent with each of the models for 
the large scale thermal structure. We also determine 
the uncertainties on the X-ray normalization and den- 
sity profile parameters due to systematic errors in the 
background and uncertainties in the cluster extent. 



which were not considered in Elbaz et al. 1995 



3.1. Temperature Profile 

Independent constraints on the IC gas thermal 
structure can be derived from the combination of the 
X-ray derived density distribution and S-Z surface 
brightness. To facilitate this, we parameterize the 
thermal structure in terms of a hybrid model for the 
radial temperature profile of the IC gas. The model 
consists of an isothermal central region extending to 
0iso beyond which the temperature decreases accord- 
ing to a polytropic model with index 7, 



Tei0) = 



n. 



l + (9/9c) 



i+{Oi,o/ear 



1/3(7-1) 



< 



> 



(23) 

This model was originally introduced by Hughes et 
al. (1988) in their study of the Coma cluster which 
was, prior to ASCA, the only cluster for which the ra- 
dial temperature profile had been measured at large 
distances from the cluster center. In Coma, the IC 
gas temperature decreases significantly at large radii. 
In the analysis of EXOSAT, Tenma, and GINGA 
data for the Coma cluster, Hughes et al. (1993) set 
7 = 1.55 and determined a best fit isothermal an- 
gular radius 9iso ~ 3.5 9c. We have generalized this 
model by allowing both 7 and 9iso to vary. 

This parametric model provides a convenient an- 
alytical description of the temperature profile mea- 
sured by ASCA. In Fig. S, we plot the approximate 



70 confidence interval for the hybrid model param- 
eters (7, 9iso) allowed by the ASCA data. These pa- 
rameters are constrained by computing the emission- 
weighted average temperature in the three ASCA ra- 
dial bins for each of the hybrid models considered. 
We determine the x^ foi' the fit of the models to the 
measured ASCA/GIS-I-SIS temperatures, leaving the 
central temperature as a free parameter. A range of 
model parameters for the thermal structure provide 
a good fit to the ASCA data; we adopt 6iso = 4.0 9c 
and 7 = 2.0 as representative values. 
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Fig. 3.— Fits of the ASCA and S-Z data to the hy- 
brid thermal models. The hatched regions contain the 
allowed values of the model parameters at 68% confi- 
dence. The ROSAT fits constrain 7 < 2.5, especially 
for small 9iso ■ The "X" marks the adopted hybrid 
model. 

3.2. Central and Average Temperature 

The discrepancy between the ASCA and GINGA 
determined temperatures makes the comparison of 
the isothermal and non-isothermal treatments diffi- 
cult. In this section, we attempt to remedy this situ- 
ation by determining a consistent X-ray temperature 
from the two experiments. 

Markevitch et al. (1996) simultaneously fit the 
ASCA/SIS-fGIS and GINGA spectra to determine 
13.3, 13.3, and 3.8 keV in the three radial bins de- 
fined in Section |3j. For the analysis including the 
thermal structure, we adopt Tco = 13.3 keV as the 
central IC gas temperature. We determine the uncer- 
tainty of the central temperature taking into account 
systematic difference between the two experiments. 
Crudely scafing the error bars of the ASCA/GIS-fSIS 



central temperature to 68% confidence, we find Teo « 
12.2^^7^ keV. Assuming the thermal structure deter- 
mined by ASCA, the central temperature is ~ 7% 
higher than the isothermal value. In this case, the 
central temperature determined by GINGA increases 
to Teo ~ 15.6l;gkeV. We adopt error bars for the 
combined GINGA and ASCA temperatures that en- 
close the 68% confidence intervals from each of the 
experiments considered individually: ASCA {T^q > 
12.2 - 0.7) and GINGA (Teo < 15.6 + 0.6). There- 
fore, assuming the ASCA determined thermal struc- 
ture, we find Teo « 13.3l^;gkeV at 68% confidence. 
In Section |7j, we use this value of the temperature 
to determine Hq under the assumption of the ASCA 
thermal structure. 

We also determine a consistent temperature for the 
case in which the IC gas is assumed to be isother- 
mal. Using the ASCA+GINGA central tempera- 
ture and the the ASCA/GIS+SIS thermal gradient, 
we estimate the emission-weighted average temper- 
ature to be Te « 12.4 keV. The ASCA/GIS-fSIS 
emission- weighted temperature is Te ~ 11.2±0.7 keV 
at 68% confidence. We adopt error bars for the 
combined GINGA and ASCA analysis that enclose 
the 68% confidence intervals from both the ASCA 
(te > 11.2-.7keV) and GINGA (te < 14.6-1-0.6 keV) 
isothermal spectral analysis. Therefore, assuming the 
IC gas to be isothermal, we find tg « 12.4t^g keV at 
68% confidence. In Section |6], we use this value of the 
temperature to determine Hq assuming the IC gas to 
be isothermal. 



3.2.1. X-ray Emissivity 

The temperature uncertainty contributes to the 
uncertainty of the central specific emissivity (Aeo) in 
equation El^. Aeo is obtained by convolving A{Te,E) 
with the ROSAT/PSPC response (|Elbaz et al. 19951). 



The emissivity varies by AAeo/K^ —~to % for 
Teo = 13.3l^'g keV; similar errors are obtained for 
an isothermal plasma with te — 12.4^^gkeV. 

We have fixed iVn and th e metallicity to th eir best 
fit value s as determined by Elbaz et al. 1995 (as was 
done by Markevitch et al. 1996| in their ASCA anal- 
ysis). However, there are uncertainties in these val- 
ues; from the ROSAT/GINGA analysis, log(7VH) = 
21.22 ± 0.03 and the metalhcity is 0.4 ± 0.1 at 68% 
confidence. The corresponding uncertainties on Aeo 
are ±2.5% and ±1%, respectively. We have assumed 
that these additional uncertainties are independent 



from the uncertainty in the temperature. In particu- 
lar, the A^H value is constrained by the PSPC spec- 
trum which is insensitive to the exact temperature. 

Adding the uncertainties in quadrature, we obtain 
a statistical uncertainty in the central specific emis- 
sivity, AAeo/Aeo =^4 8 %• Tliis Uncertainty is corre- 
lated with the temperature uncertainty (Aeo decreases 
with increasing temperature) and is included in the 
uncertainty in Hq due to variations in temperature. 

3.3. X-ray Surface Brightness Analysis 

3.3.1. Isothermal Model 

The normalizations {Nx, Nsz) and the density 
model parameters (/3, 9c) are correlated. The term, 
Ng^/Nx, which appears in equation |2l|, must be de- 
termined consistently from the X-ray and S-Z data. 
The statistical precision and spatial resolution of the 
X-ray imaging data for A2163 are much better than 
the corresponding precision in the S-Z data. As a re- 
sult, the density profile shape is constrained by the X- 
ray surface brightness profile. The error on the peak 
Comptonization yo, determined from the S-Z data, is 
dominated by the statistical errors on the S-Z data 
and not by the uncertainties in the density model pa- 
rameters. We therefore express Ng^/Nx in terms 
of the nearly independent quantities N^^ (determined 
from X-ray imaging data only) and yo, 
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(24) 



(25) 



bxo is the central X-ray surface brightness, and B is 
the incomplete beta function. 

Because bxo and the density model parameters are 
not independent, we directly compute the uncertainty 
in N^ from fits to the original PSPC surface bright- 
ness profile obtained by Elbaz et al. (1995). A x^ 
procedure, as used in Avni (1976), is applied to de- 
termine the uncertainty with the other parameters 
(/?, 9c, and background level) being optimized for any 
given value of N-^. We find the statistical uncer- 
tainty to be AN'x/N'x = ±3.5% at 68% confidence. 
The absolute calibration of the PSPC is uncertain 
by ~ ±5%. Adding the uncertainties in quadrature, 
we determine the X-ray normalization to be uncer- 
tain by AN'x/n'x = ±6.1%. At 68% confidence. 



the density model parameters are (3 — .eiel^Jg and 
6'c = 1.20 ±0.05'. 

The derived density profile shape depends on the 
background level. The background level is a free pa- 
rameter in the fit to the surface brightness profile 
and the statistical errors in the background are in- 
cluded in the uncertainties of the density model pa- 
rameters. In the above analysis, variations of the 
background with position due to improper vignetting 
corrections and contributions from unresolved point 
sources (especially at the outer part of the field of 
view) are not taken into account. As discussed in El- 
baz et al. 1995, these background variations are ^ 9% 



Squires et al. (1996) determined the effect of the back- 
ground uncertainty on the determination of the den- 
sity model by repeating the fits with the background 
±10% from the nominal value. The additional un- 
certainties in the model parameters are A/3 = ±0.03 
and A6'c — ±0.1'. The corresponding additional un- 
certainty on the structural form factor (©520/6x0) is 
±2.5%. Adding these systematic errors in quadrature 
with the statistical errors, we obtain [3 = .616 ±0.031, 



Oc = 1.20 ± .11', and AN^/N^ 



-x 



±6.6%. 



3.3.2. Truncated Model 



The isothermal (i model (equation |2^) is a good 
fit to the X-ray data within the radius of maximum 
detection {Omax = 15 9c)^ at which the s ignal to noise 
ratio drops to unity ( Elbaz et al. 1995 ). For values 
of /3 < 1.0, this simple model implies an infinite mass 
of IC gas. Therefore, at large radius, the density dis- 
tribution of A2163 (/3 = .616) must decrease more 
rapidly than the inferred density profile. One way 
to obtain a more realistic distribution is to assume 
that the gas density falls to zero at a cut-off radius 
[Ocut] where 9 cut > 9 max- We have fit the X-ray data 
with models for the surface brightness due to a trun- 
cated density distribution for a range of cut-off radii, 

"max — ^cut -1 C50. 

In the most extreme case {9cut = 9 max), we find 
(3 = .625 and 9c = 1.22'. These values are both within 
the 68% confidence intervals determined using the 
standard j3 model. The corresponding change in the 
X-ray normalization is Nx{9cut — 15 0c)/A^x(oo) = 
.98. 

3.3.3. Hybrid Model 

In the presence of thermal structure, the observed 
morphology of the X-ray surface brightness will be a 



function of both the temperature of the IC gas and 
the energy band over which the measurement is made. 
We repeat the analysis of the ROSAT/PSPC data 
assuming models for the thermal structure given by 
equation p3[ 

We first compute the specific emissivity by con- 
volving K{Tc, E) with the ROSAT/PSPC response in 
the considered energy band. A redshifted isothermal 
plasma model is used, where the hydrogen column 
density and the heavy metal abundance are fixed to 
their best fit isothermal values. Assuming values for 
the density and thermal model parameters, we then 
create a model for the X-ray surface brightness by 
numerically integrating equation O. Each model, af- 
ter convolution with the PSPC PSF and the addition 
of a constant background, is fit to the observed sur- 
face brightness profile. For a given thermal model, 
we then determine the best fit density model param- 
eters and X-ray normalization, by minimizing the x^ 
of the model fit. We have done this for a grid of 
hybrid thermal models described by the parameters 
{9iso, 7). For the adopted thermal model (7 = 2.0 and 
6»,3<, = 4.0 6*^), we determine /3 = 0.640 and 6*^ = 1.26'. 
These values are listed, along with the isothermal 
model parameters in Table |l|. In Section [t], each 
model for the thermal structure is combined with its 
corresponding ROSAT/PSPC derived density profile 
in order to determine a consistent model of the S-Z 
surface brightness. 

The changes in the best fit density profile param- 
eters are, for most models considered, smaller than 
the corresponding statistical uncertainties. However, 
for all models with values of 7 > 2.5, especially those 
with small 9iso, the fit to the PSPC surface bright- 
ness is bad. The temperature reaches such a low 
value that the X-ray emission in the PSPC band cuts 
off sooner than the observed surface brightness distri- 
bution. This constrains the acceptable temperature 
model parameters to lie within the range plotted in 
Fig. I 

4. S-Z Observations 
4.1. Instrument 



All S-Z observations were made using the Sunyaev- 
Zel'dovich Infrared Experiment (SuZIE) bolometer 
array at the Caltech Submillimeter Observatory 
(CSO) on Mauna Kea. A detailed description of 
t he SuZIE instrument h as been presented elsewhere 
( Holzapfcl et al. 1996a ). SuZIE is a 2 x 3 array of 



Temperature and Density Model Parameters 



Thermal Structure 


kT,o [keV] 


7 


&iso [^c] 


P 


ec 


Isothermal 


UAtli 


1.0 


- 


0.616 


1.20 


Truncated 


UAtli 


CX) 


15.0 


0.625 


1.22 


Hybrid(ASCA) 


13.3tli 


2.0 


4.0 


0.640 


1.26 



Table 1: Temperature and density model parameters for the isothermal, truncated, and hybrid models. 



300 mK bolometric detectors optimized for the obser- 
vation of the S-Z effect in distant {z > .1) clusters of 
galaxies. A tertiary mirror re-images the Cassegrain 
focus of the telescope to an array of parabolic con- 
centrators which couple the radiation to the six indi- 
vidual composite bolometers ( Alsop et al. 199^ ). In 
order to reduce spill-over, a 2 K stop limits the illu- 
mination of the CSO 10.4 m primary mirror to 50% 
of its total area. The array consists of two rows sep- 
arated by 2.2'; each row consists of three co-linear 
array elements separated by 2.3'. Each array element 
produces a beam on the sky which is approximately 
1.75' FWHM. 

The spectral response of the array elements is de- 
termined by a common set of metal-mesh filters po- 
sitioned between the parabolic concentrators and the 
cold stop. The filters can be changed for observations 
in several mm- wavelength passbands. The 2.1 mm fil- 
ter passband is designed to maximize the ratio of the 
S-Z signal to the sum of the atmosphere and detector 
noise and has a FWHM A A/ A = .11. 

Array elements within each of the two rows are 
electronically differenced by placing pairs of detec- 
tors in AC biased bridge circuits. The output of each 
bridge is synchronously demodulated to produce a 
stable DC voltage proportional to the instantaneous 
brightness difference between the two array elements. 
The electrical bias on each detector is adjusted so 
that the responsivities of the three detectors in each 
row are matched ( Glezer et al. 199^ ). The bolometer 
differences strongly reject signals common to both el- 
ements, due to fluctuations in the temperature of the 
300 mK heat sink and atmospheric emission. Each 
row produces three differences, two with beams sepa- 
rated by 2.3' and one with a 4.6' separation. During 
all observations, the telescope is fixed in place and 
the rotation of the earth drifts the source across the 
array of detectors. Keeping the telescope fixed while 
taking data eliminates signals due to modulation of 



the telescope's sidelobes and microphonic response of 
the detectors. Between scans, the array is rotated 
about the optical axis so that the two rows of the ar- 
ray are kept parallel to the direction of the scan. The 
difference signal from each pair of detectors and the 
absolute voltage on each of the detectors is low-pass 
filtered at 2.25 Hz and digitally sampled at 5 Hz. 

4.2. Observations 

A2163 was first observed on the nights of April 
23-26, 1993 for a total of 16 hours of integration at 
2.1mm; these observations have been described in 
W94. W94 did not include the data from the night of 
April 23 in their analysis due to the short observation 
time and comparatively poor quality of the data. We 
include these results in this reanalysis, however, they 
make almost no contribution to the weighted mean. 
We observed A2163 at 2.1mm on the nights of 1994 
April 4, 9, 10 and 11 for a total of 8 hours. The ar- 
ray was positioned so one row passed over the X-ray 



center a s determined by the ROSAT/PSPC ( |Elbaz ei 
al. 1995), while the second row passed 2.2' in declina- 
tion to the south. Scans were 30' long and begun at a 
Right Ascension Offset (RAO) leading the source by 
either 12' or 18'. In Fig. |2| we show a schematic of 
the 1994 scan strategy overlaid on the measured X- 
ray surface brightness contours. First, the telescope 
is pointed to a position leading the source by one of 
the two fixed RAOs. The telescope tracks this posi- 
tion while the data from the last scan is compressed 
and stored, the reciever is rotated, and the computers 
perform housekeeping chores. After about 10 s, the 
telescope stops tracking and a new scan begins; while 
gathering data, the telescope is held fixed. At the 
end of the scan, the telescope is reset to the alternate 
RAO and the cycle is repeated. As is described in 
Section 5.3.4., alternating the RAO between sequen- 



tial scans provides a sensitive test for instrumental 
baselines. 



Over the course of the observations, the zenith an- 



gle ofA^li33rajigcdfrom50^ 



depth of the atmosphere at the CSO was simuhane- 
ously measured with a tipping radiometer operating 
at 225 GHz. Converted to the frequency and zenith 
angle of these measurements, the atmospheric absorp- 
tion ranged from 2 — 4% over the course of the obser- 
vations and has been corrected for in the calibration 
of the data. 

As a test for a possible instrumental baseline 
we observed regions free of known sources, se- 
lected so they could be observed over a similar 
range of hour angle and zenith angle as A2163. 
The 1993 observations were coupled with 6.7 hours 
of integration on a patch of sky centered at 
13'' 11" 29", -1°20'11" (J2000). In 1994, we ac- 
cumulated 8.5 hours of data on patch #1 centered 
at 10*^24" 25", 3° 49' 12" (J2000), and 7.5 hours on 
patch #2 centered at 16^32'° 44", 5° 49' 41" (J2000). 

4.3. Calibration 

In April 1993, scans of Jupiter, Saturn, Mars, and 
Uranus were made in order to map the beam shapes of 
the instrument and calibrate its responsivity. In W94, 
Mars was used as the primary calibration source. The 
scans of Mars were made with the receiver rotated 
by an angle somewhat larger than that used in the 
source observations. We have reanalyzed the April 
1993 observations using Uranus for the calibration. 
The scans across Uranus were made with a receiver 
rotation angle of 10° which this lies near the center of 
the range of rotation angles under which A2163 was 
observed. This reduces the uncertainty in calibration 
due to changes in beam shape with receiver rotation. 

In April 1994, both Uranus and Jupiter were avail- 
able as calibration sources. As in the reanalysis of 
the 1993 data, Uranus was used both to map the 
beams and as the primary responsivity calibrator. 
The beam-shapes measured with Jupiter and Uranus 
exhibited a small and reproducible dependence on the 
rotation angle of the instrument about the optical 
axis. Over the course of the A2163 observations, the 
rotation angle varied from —45° to +37°. The cal- 
ibrated models computed for the range of rotation 
angles differ by < 5%. 

In order to calibrate the data, we determine the 
ratio of the S-Z brightness to the brightness of the 
calibrator. The brightness of Uranus is found from 
a third order polynomial fit to the measured bright- 



ness temperature as a function of wavelength ( prif 
fin fc Orton 199^ ). They assign an uncertainty of 



±6% to the brightness of Uranus, most of which arises 
from the ±5% uncertainty in the absolute brightness 
of Mars from which Uranus is calibrated. The spec- 
tral response of the array elements, including detailed 
checks for out of band leaks, have been measured us- 
ing a Fourier transform spectrometer. The central fre- 
quencies of the six array elements are determined with 
an absolute accuracy of better than 1% and exhibit a 
r eproducible scatter of less than 1% about the mean 
( Holzapfel et al. 1996a ). Including uncertainties in 
the spectral calibration, measured beam-shapes and 
absolute brightness of Uranus, the uncertainty in the 
absolute calibration of the instrument responsivity to 
the thermal S-Z effect is estimated to be ±8%. 



5. 2.1mm Analysis 
5.1. S-Z Data Set 

Each of the two rows consists of three detectors: 
si, 32, and S3. Detector signals are differenced in 
pairs (in hardware) to form three difference signals: 
di2 = si - S2, ^23 = S2 - S3, and ^31 = S3 - Si. 
These differences correspond to angular chops of 2.3', 
2.3', and 4.6' respectively. Each detector contributes 
to two of the difference signals, therefore, the three 
detector differences are not independent. To remove 
the degeneracy, the two 2.3' differences are differenced 
to form a triple beam chop. 



123 



= di2 - d23 ^ Si - 2S2 + S3 



(26) 



This combination of detector signals has the added 
benefit of being insensitive to linear changes in bright- 
ness across the array, and therefore provides supe- 
rior sensitivity in the presence of varying atmospheric 
emission. The data set consists of 4 difference signals: 
c^3ii ti23, ^64 J a-iid ^456 Corresponding to the 4.6' dif- 
ference and triple beam chop for each of the two rows 
of detectors. From now on, we refer to these four 
difference signals as dk, where k ranges from 1 to 4. 
The average signal of the undifferenced detectors is 
also computed, 






Sk 



(27) 



We use the average single detector signal as a monitor 
of the absolute atmospheric emission and to correct 
any residual common-mode response of the detector 
differences. 
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We clean the raw data of transients due to the 
interaction of cosmic rays with the detectors. The 
bolometer and electronic time constants are fast 
enough that the system recovers within < 1 s. These 
transients are idetified using an algorithm that com- 
putes the derivative of the data in the scan and then 
looks for the large positive and negative slopes asso- 
ciated with a cosmic ray event. Less than 5% of the 
data are identified as contaminated by cosmic rays. 
The raw data are then binned into 3 s bins corre- 
sponding to 15 5 Hz samples or 0.75' cosJ on the sky. 
Samples flagged as bad due to cosmic rays are left out 
when the bin averages are computed; bins with less 
than 8 samples are not included in the analysis. 

5.2. S-Z Surface Brightness Model 

The computed S-Z surface brightness profile of 
A2163 is ^ 5' FWHM, greater than the largest beam- 
throw of SuZIE, which is 4.6'. In order to accu- 
rately determine the central surface brightness, we 
must simulate the observation of an extended source. 

A model for the surface brightness of the S-Z 
thermal component is computed from the X-ray sur- 
face brightness determined density profile and the as- 
sumed thermal structure. 



IsZiy{0, (p) QC / fnh{x,Te)d^ 



(28) 



We express the surface brightness morphology in 
terms of a dimensionless form factor, 



s.. 



' /sz..(o,o) //„(o,o,e)/*.(o,o,Orfe' 

(29) 
where Si (in order to take relativistic corrections into 
account) is evaluated at the S-Z intensity weighted 
band center for each array element. The surface 
brightness model (5*^) is convolved with the single de- 
tector beam-maps constructed from scans over plan- 
ets (Vpi) to generate a model for the response of the 
detector signals. 



smK,Ti{0) = 
Alszi 1 
Ipi Up 



Vp,{0 



(30) 

\Sk.t{o', (t)')de'd4)' , 



where Vlp is the solid angle subtended by the plane- 
tary calibrator. The ratio of the relativistically cor- 
rect S-Z brightness to the brightness of the planet is 



Iszt ^ j ft{v)Iszvdv 
I Pi ^ Ji{i')Ipudv 



(31) 



where fi{v) is the measured spectral response of each 
array element and Ipi, is the intensity of the planetary 
calibrator. 

The models for the response of the array elements 
are differenced to create models for the response of 
the detector differences. The ~ .05' resolution differ- 
ential source model is binned identically to the scan 
data to determine the model signal for each of the 
~ .75' data bins. The source models are then desig- 
nated as mki{RA), where k is one of the four detector 
differences, i is the position (by bin number) in the 
scan, and RA is a ofset of the model from the nominal 
X-ray determined position. 

5.3. Peak Comptonization and Uncertainties: 
Isothermal Model 

5.3.1. Coadded Data Analysis 

We coadd the difference signals over many scans to 
create a high sensitivity scan of the differential surface 
brightness as a function of RA. For each difference 
signal (k) and scan (j), we clean the data by removing 
the best fit linear baseline. 



^kji 



dk 



kji ^kj 



ib 



kj 



(32) 



The value of each coadded bin i (xki) is given by 
the weighted sum of its values in each of Ns scans. 



E 



RMSf 



Xki 



E 



(33) 



RMST" 



Each bin is weighted by the residual RMS of the scan 
with the best fit (minimum RMS) linear baseline re- 
moved. 



N,, 



RMS^^. 



Ex2 



i=l 



kji 



(34) 



(iV,-l)' 

where Ni, is the number of bins in a scan. The uncer- 
tainty in the value of each bin is determined from the 
weighted dispersion of that bin for each scan about 
the mean of N^ scans. 
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Ns 
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Ns 
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(35) 



j=i 
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The coadded scans are fit with the isothermal 
source niodel determined in Section ^.2j . The best fit 
source amplitude (yo) and position (RA) are found 
by minimizing the x^ of the fit to all four difference 
signals to the appropriate differential source models, 



Nt 



fc=l 1=1 



{xki - 2/0 ruki (RA) - Uk -ibkY 



(36) 



'ki 



where the sum over k is for all the difference signals. 
This procedure accurately determines the source am- 
plitude, so long as the signal is a small contribution 
to the structure in the individual scans. If the source 
contributes significantly to the structure in a single 
scan, a scan with less signal will have lower RMS and 
will be weighted higher in the coadd. Coadding the 
scans using equations p3-p4 could introduce a bias 
in the amplitude of the coadded scans. However, the 
position determined from the coadded data should be 
accurate. Once the source position is known, we can 
repeat the analysis in a way that takes the contribu- 
tion of the source to the scan weighting into account. 
The construction of the coadded scans is repeated 
with the bins weighted by the residual RMS of the 
scans after the removal of a the best fit source model, 
residual common-mode signal, and linear baseline, 

RMS^^. = (37) 

Y. {'^kji " Vkj mki{RA) - akjSji - akj - i buj) 

i=l 

{N, - 1) 

Allowing the source model to vary, while determining 
the weights, ensures that the weighting does not bias 
the signal amplitude. In Appendix A., we discuss 



the removal of the residual common-mode signal and 
demonstrate that it has no systematic effect on the fit 
results. We then coadd the binned scans with the best 
fit linear baseline and residual common-mode signal 
removed. 



•^kji 



dk 



kji ^kj ^ji ^kj 



ik 



kj ' 



(38) 



Using equations p3^ and Bq, we recompute the bin av- 
erages and uncertainties which are free of bias due 
to the non-zero source amplitude. In Fig. 0, we show 
coadded data scans for the 1994 observations of A2163 
and a patch of blank sky. 

The determination of accurate uncertainties for the 
coadded scan fits is difficult due to the presence of 
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X-ray Center 




Blank Sky 
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-10 10 

RA offset from X-ray center [arcmin] 

Fig. 4. — Coadded 4.6' difference data for scans across 
the X-ray center of A2163 , 2.2' to the south, and re- 
gion #2 of the 1994 "blank sky" . The data is plotted 
in terms of a Rayleigh- Jeans temperature difference 
completely filling one of the beams. The heavy curve 
is the best fit isothermal model to the 4.6' and TBC 
data, and the light curves correspond to the ±lcr un- 
certainties on the model amplitude. 

atmospheric emission which results in noise corre- 
lated between bins and detector differences. In Sec- 
we use the fact that the noise from one 



5.3.2 



tion 

scan to the next is uncorrelated to determine the un- 
certainty in the total data set and provide an accurate 
calibration of the uncertainty in the coadded fits. 

5.3.2. Single Scan Fits 

Once the source position is known, each scan can 
be used as an independent measurement of the source 
amplitude. The source amplitude for each scan is 
found by minimizing the x^ of the model fit to the 
desired data scans. 



k i=l 



(39) 



yoj mkijRA) ~ UkjSji - akj - i bkjY 



where RMSfcj is given by equation pSl The uncer- 
tainty in the peak Comptonization (cfyj) is simply 
the change in tjqj corresponding to Ax^ — 1. We 



12 



determine the mean peak Comptonization for a given 
observation by averaging the values j/oj for each of the 
Ns scans weighted by ayj , 



2/0 



Ns 



(40) 



The weighted dispersion of the scan amphtudcs about 
the mean is used to estimate the uncertainty in the 
determination of the mean for each observing period, 





\ 


Ns 1 .2 


CTy 





(41) 



In Fig. 1^, we show a histogram of the peak Comp- 
tonization results for simultaneous fits to both array 
rows for every A2163 scan taken in 1993 and 1994. 
The mean and width of this distribution can be used 
to determine the peak Comptonization and uncer- 
tainty in the simple case when all scans are weighted 
equally. Using this method, the mean is unchanged 
while the uncertainty in yo is ^ 30% higher than that 
found when the scan amplitudes are weighted using 
equation n^. We do not use the results of this simple 
analysis, however, it illustrates the independence of 
the results from the details of the analysis. 

For each night, RAO, and scan declination, we 
use equations p9|-pl| to compute the mean amplitude 
and uncertainty; these results are listed in Table H. 
The single scan amplitudes are combined to deter- 
mine an average peak Comptonization for all the data 
at each of the three scanned declinations; these re- 
sults are listed at the bottom of the columns in Ta- 
ble 0. Due to correlated sky noise, the results of the 
fits to each of the two rows are not independent. In 
order to estimate the uncertainty of the total data 
set, we simultaneously fit to the signals from both 
rows. Assuming the sky noise between scans to be 
uncorrelated, we use equations ^ and ^ to deter- 
mine yo = 3.73 ± .35 x 10"^. 

The average source amplitudes determined by the 
single scan fits agree with the results of the fits to the 
coadded data, while the uncertainties are a factor of 
1.4 to 2.0 higher than those determined from the fits 
to the coadded data. This is because the atmospheric 




-o.ooa 0.002 

Single Scan CompLonization Value 



Fig. 5. — Histogram of peak Comptonization values 
from simultaneous fits to all four differential signals 
in each individual scan. The distribution is well ap- 
proximated by the superimposed Gaussian, the mean 
is marked with a vertical line. 



noise, which dominates the noise in the scans, is cor- 
related between difference signals and bins in a given 
scan. In order to use the coadded data to determine 
confidence intervals for the position and morphology 
of the S-Z surface brightness, we must take the corre- 
lation of the noise into account. To do this, we scale 
the uncertainties of the coadded bins so that the un- 
certainty in the coadded scan amplitude is the same 
as that from determined the single scan fits. The 
scaling factor is insensitive to the details of the clus- 
ter position and morphology. We can then use the 
coadded data to determine the confidence intervals 
for parameters such as the model position which can- 
not be determined from the individual scan fits. All 
fits to coadded data discussed in this paper are for 
data with bin uncertainties scaled to compensate for 
correlated noise. 

5.3.3. Peak S-Z Surface Brightness Position 

With the coadded scan bin uncertainties appropri- 
ately scaled, we use a maximum likelihood indicator 
to determine confidence intervals in the model ampli- 
tude and position. We determine the position of the 
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Peak Comptonization xlO"* in A2163: Isothermal Model 



Date 


RAO 


A(5 = 0' 


AS ^ +2' 10" 


AS = -2' 10" 






Peak ( 


I^omptonization, 


yo X 10^ 


1993 April 23 


16.5' 


0.72 ±2.78 


-3.34 ±4.64 






22.5' 


2.52 ±2.15 


9.33 ±6.36 




1993 April 24 


16.5' 


5.06 ±1.39 


7.33 ±2.36 






22.5' 


3.40 ±1.66 


0.51 ±2.32 




1993 April 25 


8.5' 


4.46 ±1.17 


1.74 ±2.26 






14.5' 


3.30 ±1.22 


4.96 ±2.20 




1993 April 26 


8.5' 


4.39 ±0.82 




1.32 ±2.32 




14.5' 


6.71 ±1.08 










2.03 ±2.48 


1994 April 4 


12' 

18' 


5.38 ±2.42 
7.01 ±1.20 










11.43 ±7.33 
15.11 ±5.59 


1994 April 9 


12' 
18' 


3.20 ±0.77 
2.41 ±0.83 










3.58 ±2.20 
1.67 ±2.61 


1994 April 10 


12' 
18' 


3.08 ±1.65 
2.25 ±1.08 










-2.18 ±3.53 
3.17±3.17 


1994 April 11 


12' 
18' 


4.36 ±1.04 
3.54 ±1.23 










1.05 ±3.10 
2.03 ±3.42 


Weighted mean 




3.87 ±0.35 


3.50 ±1.10 


2.77 ±0.97 



Table 2: Results of fits of the isothermal model to the data for each observation, S, and RA. The total results for 
each declination (bottoms of columns) are calculated in the same way as the individual sets using the distribution 
of the scan amplitudes. 



peak S-Z surface brightness to be offset from the peak 
X-ray surface brightness by ARA = +.35 ± .14' at 
68% confidence. The value of the peak Comptoniza- 
tion decreases by less than 1.0% from its value at the 
central position over the allowed range of RA. 

We estimate the uncertainty in pointing the SuZIE 
instrument to be < 10 . ROSAT astrometry is typ- 
ically uncertain by 10 — 15", unless special care has 
been taken to locally align X-ray and optical frames. 
Because no X-ray bright optical sources were present 
in the ROSAT field, this was not possible for A2163. 
Including the uncertainty in the ROSAT and SuZIE 
pointing, the offset between the X-ray and S-Z peak 
surface brightnesses is then ARA = +.35 ± .33' at 
68% confidence. 

We do not have sufficient coverage to accurately 
locate the peak S-Z surface brightness in declination. 
Assuming the peak S-Z and X-ray surface bright- 
nesses to be coincident, the S of the peak surface 
brightness is uncertain by the quadrature sum of the 
ROSAT and SuZIE pointing uncertainties, AS ss .30'. 
If the true position of the peak surface brightness is 
offset from the declination of our scan, the measured 



peak Comptonization will be smaller than the true 
value. Fitting the data to a model for scans offset 
in S from the peak surface brightness, we find that 
the peak Comptonization is ~ 2.0% smaller for AS = 
±.30'. The uncertainty in the position of the peak S-Z 
surface brightness, therefore, results in an uncertainty 
in the peak Comptonization, Ayo/yo — +2.2%. This 
contribution to the total uncertainty is listed under 
"Position" in Table |[ 

5.3.4- Baseline 

We analyze the "blank sky" scans across regio ns o f 
sky free of known sources (described in Section 4^ ) 
exactly as the source data, to test for an instrumen- 
tal baseline. An instrumental baseline with the same 
temporal features as the source model would be de- 
tected as a non-zero Comptonization for these regions. 

In the April 1993 observations, the RAO of se- 
quential scans was alternated in the same way as the 
for the source data. In W94, this data was used 
to estimate the baseline contribution to the uncer- 
tainty in the source amplitude. Reanalyzing this 
data using the ROSAT density model and the resid- 
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ual common-niode signal subtracted scans, we find 
yo = -0.1 ±4.2 X 10-5 

In the April 1994 observations (patch #1 and patch 
#2), the RAO was alternated between 12' and 18' in 
sequential scans. The results of these observations 
have been used to place upper limits on arcminute 
scale CMB anisotropies (Church et al. 1996). An- 
alyzing the baseline data in the same way as the 
source data, we find yo = 1-3 ± 3.1 x 10"^ and 
yo = 1.8±3.1xl0^5for the two patches, respectively. 
In Fig. H, we show the coadded 4.6' difference for the 
scans on patch #2. Combining the results from all 
three patches, we determine y^ — 1.45 ± 1.9 x 10^^, 
indicating no significant instrumental baseline. Com- 
bining the results from multiple patches reduces the 
signal from potential sources of confusion but should 
accurately determine the effect of any common in- 
strumental baseline. One problem with this method 
is that the baseline data and source data, although 
gathered at similar azimuth and zenith angle, are not 
gathered simultaneously or even on the same evening. 
Any baseline signal, because it is not significant in a 
single scan, must be correlated in time between sev- 
eral scans. This does not guarantee the baseline will 
be constant over the course of an observation. 

We have devised a method which allows us to use 
the source data to test for the presence of a scan- 
correlated baseline. Differencing scans taken adjacent 
in time at two different RAOs allows the subtraction 
of the common baseline with little effect on the ex- 
pected signal. Each pair of scans are differenced by 
subtracting their raw 5 Hz sampled time streams. The 
data are then analyzed exactly as described in Sec- 
tions |]-|3a. 



We construct a model for the differential surface 
brightness by differencing two models with the ap- 
propriate RAO difference between them. We choose 
the RAO difference (6 ) to be large enough so that 
the two models are well separated in the scan and 
therefore will not subtract much signal when differ- 
enced. The differenced model is then fit to the corre- 
lation corrected coadded scans in order to determine 
the source amplitude and position. Fitting to the 
isothermal model, we find yg = 3.49 ± .40 x lO""* and 
Ai?A = -H0.19 ± 0.18' at 68% confidence. These re- 
sults agree, within the quoted statistical errors, with 
those determined from the undifferenced data. 

The lack of any significant change when the scans 
are differenced, as well as the null result for the "blank 
sky" fields, indicates that there is no significant in- 



strumental baseline with the same temporal depen- 
dence as the source model. The difference between 
the results of the fits to both RAOs and the differ- 
enced scans is treated as an additional contribution 
to the uncertainty in j/o and is listed under "Baseline" 
in Table pi Including the additional uncertainty due 
to the possibility of an instrumental baseline, the re- 



sults from Section 5.3.2. become yo 



3.73j::43 X lO-'^. 



5.4. Density Model Uncertainties 

We consider now the effect of the uncertainties in 
the density model on the determination of yo- We 
have fit the S-Z data with the range of density mod- 
els allowed by the X-ray analysis. The uncertainties 
in the density model parameters (Section 3.3.1.) con- 



tribute an uncertainty in the peak Comptonization 
of Ayo/yo =^3:6 % at 68% confidence. In W94, the 
data were fit with a template derived from the den- 
sity model determined from th e analysis of Einstei n 
data: 6'c = 1.15 and (i = 0.59 ( |Arnaud et al. 1992| ). 
In this paper, these data are reanalyzed using the 
ROSAT/PSPC density model which results in a peak 
Comptonization lower by ~ 4%. 

T he tru ncated density model, introduced in Sec- 
tion 3.3.2., is used to generate a model for the isother- 



mal S-Z surface brightness. Fitting the measured S-Z 
surface brightness with the model for the truncated 
surface brightness, we find yo{Omax)/yo{oo) = .925. 
The angular scale factor {Qsz) also decreases when 
the density model is truncated; this nearly compen- 
sates for the decrease in the measured peak Comp- 
tonization. Even in the most extreme case, the S-Z 
normalization is nearly independent of the presence 
of a density cut-off and Ns z (Smax) / Ns z {00) = .985. 
The uncertainty in the peak Comptonization due 
to the adopted density model is the sum of the uncer- 
tainties due to the uncertainty in the model parame- 
ters and the assumed functional form for the density 
profile. Adding these contributions to the uncertainty 
in quadrature, we determine Ayo/j/o =-8 3 '^^ This 
analysis is performed assuming an isothermal IC gas; 
in the presence of thermal structure of the form dis- 
cussed in Section n\ the effect of the cut-off density 
model will generally be smaller. 

5.5. Isothermal Temperature Uncertainty 

Due to our use of the relativistically correct cal- 
culation of the Comptonization (Section |2.1.| ), the ef- 
fective peak Comptonization derived from S-Z data 
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depends on the assumed temperature. For tempera- 
tures in the range allowed by the GINGA and ASCA 
spectral analysis (Section 3.2.), yo oc Tg/^ changes 



by li;2%- Including uncertainties in the baseline, cal- 



ibration, and density distribution and temperature of 
the IC gas, we detrermine yo 



3.73til X 10-4. 



6. Hubble Constant: Isothermal gas 

Using equations |2l|, ^ and |2^ we can express dA 
in terms of the nearly independent quantities N-^ , T^ 
and 2/0, 






kT, 



A, 



47r(l + z) 



^2 



(42) 



Apart from relativistic corrections, the errors in yo 
and Te can be considered as independent; the cor- 
responding contributions to the total uncertainty for 
Hq are added in quadrature. 

Assuming a standard (3 model, we determine the 
contribution of the uncertainty in the density model 
to the uncertainties in N^ (Section [3.3.1. ) and yo 
(Section 5.4.). The change in yo is small and has 
the same sign as the change in AN^; therefore, 
AHq/Hq « AN^/N-^ = ±6.6% is an overestimate 
of the true uncertainty. This contribution to the 
uncertainty in Hq is listed in Table O under "X-ray 
normalization". Compared to the contributions due 
to the uncertainties in the S-Z and central temper- 
ature results, the contribution of the density model 
to the total uncertainty is nearly negligible. Taking 
only the uncertainties in the S-Z and X-ray normal- 
izations into account, we determine Ho{qo = 5) = 
59.6ti3;ikms-^Mpc-\ 

In order to account for the systematic difference 
in the ASCA and GINGA temperatures, we use the 
estimated isothermal X-ray temperature determined 
in Section p.2j . The uncertainty in Hq due to Tg in- 
cludes the effect of the temperature uncertainty on 
A(Te) as determined in Section 3.2.1., In Table ^, 
we list the results for Hq, as a function of the as- 
sumed thermal structure and the included sources of 
uncertainty. Including the uncertainties in the S-Z 
and X-ray normalizations and X-ray temperature, we 



determine Hq^qq 



59.6 



+40.7 
-22.6 



kms ^Mpc"^ 



We also determine the effect of the truncated den- 



sity model described in Section 3.3.2. on the value 



of Hq. Because Nsz and Nx are insensitive to the 
details of the density model, so is the value of Hq; 



for the most extreme case {9cut — Omax), we find 
H(){d„iax) / Hq{oo) = 1.01. For the remainder of this 
paper we assume the IC gas density can be described 
by a standard /3 model with 9 cut = 00. Because the 
X-ray normalization contains the uncertainty due to 
density model parameters, the additional contribu- 
tion to Hq is simply that due to the assumption of 
the density functional form, AHq/ Hq = -1-1.0%. This 
contribution to the uncertainty in Hq is listed under 
"Density Model" in Table 0. This analysis is per- 
formed assuming an isothermal IC gas; in the pres- 
ence of thermal structure of the form discussed in Sec- 
tion [t], the uncertainty due to the assumption of a 
cut-off density model will generally be smaller. 

we determined the effect of the 



In Section 5.3.3 



uncertainty in the position of the peak S-Z surface 
brightness on the value of yo- The uncertainty in 
the position contributes an uncertainty to the Hub- 
ble constant of AHq/Hq = —4.5%; this contribution 
to the total uncertainty is listed under "Position" in 
Table |. 

7. Large-scale Thermal Structure 

In this Section, we repeat the determination of the 
Hubble constant under the assumption of a thermal 
structure of the form described by equation |2^. As 

we repeat the analy- 



was outlined in Section 3.3.3 



sis of the ROSAT surface brightness in order to red- 
erive the density model parameters (/3, 6c) for each 
set of thermal structure parameters (7, 9iso)- For the 
adopted hybrid model (7 = 2.0, 6iso = 4.0 0c), the 
reanalysis of the ROSAT data yields density param- 
eters /3 = 0.640 and 6c = 1.26'. Using equations [2q - 
|3l| we combine the models for the IC gas tempera- 
ture and density to create models for the S-Z surface 
brightness. Each model is fit to the coadded data 
scans in order to determine the peak Comptoniza- 
tion. In Fig. ^ we plot the 68% confidence inter- 
val in the {6iso, 7) plane for fits to the coadded S-Z 
data. The S-Z measurements lack sufhcient sensitiv- 
ity to detect a thermal gradient of the type deter- 
mined from recent ASCA results. The ratio of the 
peak Comptonization to that determined assuming an 
isothermal gas is plotted as a function of 6iso and 7 
in Fig. p. For the range of temperature model param- 
eters allowed by the ASCA results at 68% confidence 
(see Fig. ^, the values of the Comptonization span 
yo/{yo)iso = -82 ± .03. Using the procedure described 
in Section 5.3.2., we determine the peak Comptoniza- 
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Hubble Constant (Thermal Structure) 



Thermal Structure N 



X 



Ho [kms^^Mpc^^] 
Nsz Nx, Nsz, Ten Total 



Isothermal 
Hybrid(ASCA) 



78.4131-3 



59. 6^ 



-40.7 
-"-22. 6 

78Atlri 



59.6t^g| 



78.4 



+59.9 
-39.8 



Table 3: The Hubble constant determined for both isothermal and hybrid thermal structure. The results are calcu- 
lated for three different cases. In the first column under Hq, only uncertainties in the X-ray and S-Z normalizations 
are included. The second column includes the uncertainty in the IC gas central temperature. The third column 
includes additional uncertainties due to peculiar velocity, astrophysical confusion, clumped IC gas, and cluster 
asphericity. 



tion from single scan fits to the hybrid S-Z surface 
brightness model to be yo = 3.07 ± .29 x lO"*. The 
results of the fits at each declination are listed in Ta- 
ble |. 

Peak Comptonization & Uncertainty 



Source 




Uncertainty 
2/0 X 104 


Statistical 




3.07 ±.29 


Baseline 




-.27 


Calibration 




±.25 


Position 




+.07 


Density Model 




-I-.08 
-.25 


Thermal Gradient 


(ASCA) 


±.08 


Central Temperature 


+ .05 
-.04 






smtil 


Peculiar Velocity 




±.23 


Radio Confusion 




+ .25 
-.05 


Primary Anisotropics 


±.12 



Total 



3.07t-^;^ 



Table 5: Peak Comptonization and contributions to 
uncertainty using the adopted hybrid thermal struc- 
ture. 

We repeat this analysis for the RAO differenced 
data. The peak Comptonization from the single scan 
fits is 2/0 = 2.80 ± .32 x 10"''. We use the difference 
between this result and the peak Coniptonization for 




2.5 5.0 5.5 4.0 4.5 5.0 
Isothermal Core Radius, 6,,„/0^ 



Fig. 6. — Ratio of the peak Comptonization to the 
isothermal value (yo/(2/o)iso), plotted as a function of 
the hybrid thermal model parameters. 



the undifferenced scans to estimate the baseline un- 
certainty. Combining the result of the fits to both 
RAOs with uncertainties in the baseline, calibration, 
and density and temperature models, we determine 



2/0 = 3.07l;54 X 10 4. The density model uncertainty 



is assumed to be the same as in the isothermal case 
(Section ^. This is justified because the X-ray sur- 
face brightness profile is insensitive to the thermal 
structure for the considered range of thermal models. 
The individual contributions to the uncertainty in yo 
are listed in Table ||. 

We combine the X-ray and S-Z results to determine 
Hq as a function of the cluster thermal structure. In 
this analysis, we make one of two assumptions about 
the results of the X-ray observations. In the first case, 
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Peak Comptonization. yo x 10^ 



Thermal Structure AS = 0' 



AS 



-2' 10" AS 



-2' 10" 



Total 



Isothermal 
Hybrid(ASCA) 



3.87 ±0.35 
3.16 ±0.32 



3.50 ±1.10 
3.34 ±1.19 



2.77 ±0.97 
2.21 ±0.91 



3.73 ±0.35 
3.07 ±0.29 



Table 4: Peak Comptonization at each declination for the cases when the IC gas is either isothermal or has the 
ASCA determined thermal structure. The uncertainties include statistical uncertainty only. The numbers under 
total are calculated from simultaneous fits to both rows of the array. 



at 68% confidence. 



Hubble Constant & Uncertainty 



Source 



Uncertainty 
kms" Mpc~ 



X-ray Normalization 

S-Z Normalization 

Position 

Density Model 

Thermal Gradient (ASCA) 

Central Temperature 



±5.2 

+30.8 
-16.5 

-3.5 

+.8 

+ 1.5 
-1.2 
+43.8 
-21.9 



Peculiar Velocity 
Asphericity 
Clumping 
Radio Confusion 
Primary anisotropics 



78.4l^H 



+ 13.0 
-10.4 

±21.7 



+2.6 
-11.4 
+6.7 
-5.9 



Total 



78.4ti-| 



Table 6: Hubble Constant and contributions to the 
uncertainty using the adopted hybrid thermal struc- 
ture. The S-Z normalization includes uncertainties 
due to statistical uncertainty, baseline, and calibra- 
tion. 



we assume that the central temperature of the IC gas 
has been determined, this requires spatial resolution 
in the spectral measurement. For the allowed range 
of thermal models, we find Nx/{Nx)iso = .96 ± .02 
and Nsz/iNsz)iso = -92 ± .02. In Fig. |, we plot the 
ratio of the true value of Hq to that determined as- 
suming the cluster gas to be isothermal at the central 
temperature. The ASCA constraints on the thermal 
structure (see Fig. ^ limit Ho/{Ho),so = 1.13 ± .02 




2,0 2.5 3.0 3.5 4.0 4.5 5 
Isothermal Core Radius, 0.^,JB,. 



5.5 



Fig. 7. — Ratio of the true Hubble constant to the 
isothermal value {Ho/{HQ)iso) as a function of the 
hybrid thermal model parameters. The value of the 
central temperature is assumed to be constant. 

We also calculate Hq as a function of thermal 
structure for the more realistic case in which only 
the emission-weighted average X-ray temperature is 
known. The values of the normalization factors 
{Nx, Nsz) are nearly unchanged from the previous 
case. Assuming the ASCA thermal structure, the cen- 
tral temperature is about 7% higher than the emission 
weighted average value. This leads to a larger offset 
in Hq from the isothermal value than when the cen- 
tral temperature is held fixed. In Figure ^, we plot 
the ratio of the true Hubble constant to that calcu- 
lated assuming the cluster gas to be isothermal at the 
emission-weighted average temperature. The ASCA 
constraints on the thermal structure (see Fig. 0) limit 
Ho/{Ho)iso = l-29l;^Q at 68% confidence. Inagaki 
et al. (1995) have investigated the effect of thermal 
structure on the value of Hq determined from the S- 
Z effect. For a simulated "Coma-like" cluster, they 
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found Hq to be underestimated by > 20% when the 
IC gas was assumed to be isothermal. They identify 
uncertainty in the temperature profile as the most se- 
rious systematic error in the value of Hq and as a 
possible source of the "low" values sometimes found 
with the S-Z effect. 




2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 
Isothermal Core Radius, (5;„„/0^ 



Fig. 8. — Ratio of the true Hubble constant to the 
isothermal value (i?o/(-ffo)iso) as a function of the 
hybrid thermal model parameters. The value of the 
emission- weighted average temperature is assumed to 
be constant. 

We have determined Hq under the assumption 
of the ASCA/GIS-fSIS thermal structure modeled 
by the adopted hybrid model (Table |]). We find 
Ho{qo = ^) = 78+^^kms"^Mpc-^ at 68% confidence, 
where we include uncertainties in the normalizations 
{Nx and Nsz) and density and temperature models. 
The uncertainties in the X-ray normalization and the 
density model are considered to be the same as in the 
isothermal case and all errors are added in quadra- 
ture. The individual contributions to the uncertainty 
in Hq are listed in Tables |3| and g. 

8. Additional Uncertainties 

The uncertainties discussed so far are associated 
with the precision of the measurements required for 
the determination of the Hubble constant. We iden- 
tify several additional sources of systematic uncer- 
tainty associated with assumptions we have made 
in order to determine Hq from the measured data. 
Among these are deviations from spherical symmetry, 
peculiar velocity, IC gas clumping, and astrophysical 
confusion. The contribution of each of these addi- 



tional sources of uncertainty to tjq and Hq are listed 
in Tables ^ and ^ respectively. 

8.1. Deviations Prom a Spherical Gas Distri- 
bution 

The assumption that the extent of the IC gas along 
the line of sight is the same as is measured in the plane 
perpendicular to the line of sight is fundamental to 
this method for determining Hq. If the distribution 
of the IC gas is prolate (oblate) with its dimension 
along the line of sight longer (shorter) by a factor Z 
than the average of the dimensions perpendicular to 
the line of sight, then the derived Hq will be modified 
from the true Hq by a factor 1/Z (BHA). The aver- 
age Hubble constant derived from measurements of 
N clusters, Hq oc ^ J2n=i ;^' ^^^^ ^^ weakly biased 
towards higher Hq depending on the distribution of 
ellipticities contained in the sample. In this paper, 
we quantify the contribution of asphericity to the un- 
certainty in Hq determined from measurements of a 
single cluster. We estimate this contribution from 
the statistics of the ellipticity of the observed X-ray 
isophotes of a sample of other clusters. 

McMillian et al. (1989) determined the ellipticity of 
the X-ray isophotes for a sample of 49 clusters. The 
average ellipticity for this sample is e = .277. This 
average has been computed without the removal of 
any of the clusters containing obvious substructure. 
For a sample of clusters, the average elongation or 
shortening along the line of sight will be smaller than 
the average deviation from circularity of the projected 
X-ray isophotes. Departures from sphericity of the IC 
gas distribution, therefore, contribute an uncertainty 
of < ±27.7% to the determination of Hq from a single 
cluster. 

8.2. Peculiar Velocity 

If the IC gas has a bulk velocity (vp) with respect to 
the CMB rest frame, there is an additional kinematic 
component to the S-Z effect. 



AIk = -Iq- 



C 



In general, the ratio of the brightnesses of the kine- 
matic and thermal components of the S-Z effect in- 
tegrated over the 2.1mm band is small. For A2163 
(re = 12.4keV), 



jAlTf{v)dy 



.139 



lO'^ kms 



(44) 
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where Vr = v^- dl is the projection of the average clus- 
ter velocity onto the line of sight, and f{v) is the av- 
erage transmission of the SuZIE system in the 2.1 mm 
band. 

We have made additional measurements of A2163 
with filter bands centered at 1.4 and 1.1 mm. Com- 
bining the mm-wavelength measurements with Teo 
we are able to limit \vr\ < 1500 kms~^ for A2163 
dHolzapfel et al. 1996b|). Bahcall et al. (1994) make 



theoretical predictions of the one-dimensional RMS 
cluster peculiar velocities for a variety of cosmolog- 
ical models. The results range from 268 kms^ for 
f2 — 0.3 cold dark matter to 614 kms^^ for 51 = 1.0 
hot dark matter. We expect the S-Z results to reach 
this level of accuracy soon; until they do, we adopt 
these theoretical limits to constrain the peculiar ve- 
locity. For the largest predicted peculiar velocities, 
we find Ayo/yo = ±7.4% and AHa/Ha ^tlH %■ 

8.3. Isobaric Inhomogeneities 

Clumping of the IC gas on sufficiently short scales 
could escape detection with the ROSAT/PSPC. Hy- 
drodynamical simulations of subclustering by Inagaki 
et al. (1995) indicate that clumping could result in an 
overestimation of Hq by as much as 10 — 20%. Here we 
consider to what extent the existence of such clumps 
could effect the measurement of Ho for A2163. 

We assume a simple model for the clumping con- 
sisting of a two-phase medium in pressure equilib- 
rium: a hot phase of temperature kTh, and cooler 
clumps overdense by a factor B with respect to the 
hot phase. The clumps are assumed to be uniformly 
distributed in volume with filling factor /. In princi- 
ple, such clumping can be constrained by X-ray spec- 
tral data; this is equivalent to placing limits on the 
superposition of two isothermal models with the same 
form for their average density profile and tempera- 
tures kTfi and kTc — kTh/ B. We define a normaliza- 
tion factor Nc oc n^p fcT"^, where rieo^ is the central 
electron density of the hot phase. Nc is normalized to 
the best fit value obtained with the isothermal model; 
in the presence of clumping, Hq/Hqjso — Nc- 

We fit the GINGA -I- ROSAT/PSPC spectral data 
with this model, leaving /, B, kTh, and Nc free to 
vary. The hydrogen column density and the iron 
abundance are also left free to vary, but are con- 
strained to have the same value for the two phases. 
The relative normalization between the PSPC and 
GINGA data is fixed to the best value obtained from 



the isothermal model. The best fit is obtained for a 
nearly isothermal model: kTh — 15.5 keV, B — 1.15, 
and Nc = 0.98, with no significant improvement in 
the fit as compared to the isothermal model. The 68% 
confidence contours for the two parameters kTh and 
B are plotted in Fig. H, with lines corresponding to 
equal / values superimposed. In Fig. [I^ we have plot 
the 68% confidence contours for the two parameters 
kTh and Nc] it is apparent that the data prohibit low 
values of Nc- 



O 



m 

o 
o 



-I — I 
d 
-a 

CD 
> 

o 



(J) 
00 
l\ 

CD 

IT) 







' 1 1 ' 

1 0.001 

1 






- 




1 

' 0.003 
1 

1 
/ 




- 






/ 
■ 
/ 

/■ 0.01 






(D 




/ 

/ 

( 

1 : • "^ 


0.03 


o.icr 


O 


1 
1 


/ 


/-<^ 


0.30" 




ii\ y ^ 


^^^^ - 


0.99 


— 


\\\ , x^ ^^ 


'^ ^>^^^ 


— 




! \ Z^^-'"'''^ ^ ' 


^^^i-^"^^ 






fC^"^^^! I 1 


sotherma 
1 1 1 





10 20 30 40 50 
Temperature kTj^ (keV) 

Fig. 9. — Contour map for the two parameters kTh 
and B of the isobaric clumping model. kTh is the 
temperature of the hot phase and B is the overden- 
sity factor of the clumps. The heavy line represents 
the limits at 68% confidence, and is determined by 
fitting the GINGA + ROSAT/PSPC spectral data. 
Contours of equal filling factor / are indicated by the 
broken lines. The best fit point is indicated by a cross. 

In the extreme case of large kT/j values, the gas 
must contain clumps at temperature ~ kTiso embed- 
ded in a significantly hotter phase. In the GINGA 
energy band, the eniissivity is a weak function of the 
temperature for kT > 20 keV; therefore, it is difficult 
to constrain high values for kTh- However, large val- 
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Fig. 10. — Contour map for the two parameters kTh 
and Nc- The normahzation factor, Nc^ is defined 
such as Nc = HqJHq/iso- A heavy line marks the 
Hmits at 68% confidence. The best fit point is indi- 
cated by a cross. 



ues of kTh are unhkely. Unless electron thermal con- 
duction is strongly suppressed by IC magnetic fields, 
the clumps would evaporate on a time-scale much 
shorter than the Hubble time. The evaporation of 
clumps in a hot phase has been studied by Cowie & 
McKee (1977) and Balbus & McKee (1982). When 
the scale of the thermal structure is less than the 
electron mean free path, the evaporation is governed 
by the saturation parameter, ctq. At a redshift of 
z = .2 and H^ = 50kms~^Mpc~^ one half of the 
ROSAT/PSPC PSF FHWM is 50kpc. For clump 
radii smaller than this scale and density less than 
6 X 10~^cm~'', the conduction is saturated (o-q > !)■ 
The evaporation time is then given by 



Re 



^evap 



='■'>"''' ^{^) '''(^M^r°^""^^^' 



\ -0.209 
^0 



(45) 
which is always much shorter than the Hubble time. 

Considering only the effect of clumps colder than 
the bulk of the IC gas, we determine A_ffo/-f^o,iso = 



— 10.0% at 68% confidence. We have not treated the 
case in which the modeled IC gas also exhibits large 
scale thermal structure; we expect the restrictions on 
the presence of cool clumps to be even tighter in this 
case. It must be pointed out that these limits ap- 
ply only to IC gas which is static and free of strong 
magnetic fields, IC gas that does not satisfy these 
conditions could be considerably more clumped. 

8.4. Astrophysical Confusion 

Although astrophysical confusion from randomly 
distributed sources is expected to be small at mil- 
limeter wavelengths ( Fischer fc Lange 1993| ), the pos- 
sibility of confusion contributing a considerable sys- 
tematic error to the determination of j/q must be 
considered. In this region of the sky, the ampli- 
tude of anisotropic interstellar dust emission calcu- 
lated from the IRAS 100 /Ltm map and scal ed with 
the sky-ave rage spectrum of dust emission ( Wright 
et al. 1991 ) is negligible. A VLA search towards 
A2163 shows evidence of a radio source 0.8' west of 



the cluster center with an inverted spectrum (Herbig 
fc Birkinshaw 1992 ). For this source, the flux rises 
from 1 to 3mJy between 6 and 2 cm, suggesting a 
flux as large as 30mJy at 2.1mm. However, Fischer 
& Radford report an upper limit of 5 mJy {2a in a 
20" 5, 10" RA beam) on point source emission within 
1' of the X-ray center at a wavelength of 3.3 mm. For 
the adopted hybrid model, a 5 mJy point source at 
2.1 mm would cause the underestimation of the peak 
Comptonization by Ayo = -1-0.24 x 10^^. The sign re- 
flects the fact that the only known radio source in the 
scan with a flux and spectrum such that it could con- 
tribute significantly is sufficiently close to the S-Z cen- 
ter that it would decrease the size of the measured S- 
Z decrement. A2163 exhibits the brightest radio halo 
yet discovered (Herbig & Birkinshaw 1995). From 



measurements at 1.5 and 4.9 GHz, the integrated flux 
from the radio halo is estim ated to be les s than 1 mJy 
in the SuZIE 2.1 mm band ( [Herbig 1995|) . This could 
contribute at most Aj/o = ±.05 x 10~^ to the peak 
Comptonization. We can use these results to place a 
conservative limit on the contribution of radio confu- 
sion to the uncertainty in the peak Comptonization, 
Ayo<g^xlO-4. 

It is also possible for the measurement of the S-Z 
effect to be confused by the presence of primary CMB 
anisotropics. Haehnelt and Tegmark (1996) estimated 
the confusion limits from primary anisotropics to 
the determination of cluster peculiar velocities. We 
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use their results to determine the effect of primary 
anisotropies on the measurement of the peak Comp- 
tonization. For f2 = 1 (CDM) models with flbaryon = 
.01 — .1, cluster optical depth r = 0.015, and the 
SuZIE beam size, they find jAtiped < 300kms~ . 
From this, we determine that primary anisotropies 
contribute an uncertainty of Ayo/j/o = ±3.6% to the 
peak Comptonization parameter in A2163. 

9. Conclusion 

We have confirmed the previous detection of the 



S-Z effect in A2163 at 2.1 mm (Wilbanks et al. 1 



The 1993 and 1994 data sets have been analyzed using 
a relativistically correct treatment for the S-Z effect 
and the ROSAT derived density profile. Assuming an 
isothermal IC gas at the estimated GINGA+ASCA 



isothermal temperature, T^ = 12.4l^gkeV, and in- 
cluding uncertainties in the fit amplitude, baseline, 
calibration, position, and density model, we find j/o = 
2>.litfi X 10-''. Combining the S-Z and X-ray data 
and including uncertainties due to the S-Z and X-ray 
normalizations and isothermal IC gas temperature, 
we find Hn{qo = \) = 60^23 kms'^Mpc^^ 

There are indications that the IC gas in A2163 is 
not isothermal. Recent ASCA results suggest a dra- 
matic decrease in the temperature of the IC gas with 
radius (Markevitch et al. 1996). The model used by 
Hughes et al. (1988) to fit the thermal structure in 
the Coma cluster is generalized and used as a two pa- 
rameter fit to the temperature profile in A2163. We 
use the combined analysis of the S-Z and X-ray sur- 
face brightnesses to place limits on models for the 
thermal structure of the gas. This is the first appli- 
cation of, what is potentially, a powerful probe of the 
state of the IC gas. So far, the S-Z measurements 
lack sufficient sensitivity to distinguish between the 
thermal structure indicated by the ASCA results and 
an isothermal IC gas. 

Adopting the ASCA thermal structure, a joint 
analysis of the GINGA and ASCA/SIS-hGIS ther- 
mal structure yields a central temperature, Teo = 
13.3^2'7keV. Using this thermal structure, and in- 
cluding uncertainties due to the baseline, calibration, 
and density and temperature models, we determine 
yo = 3.07t:54 X 10--*. Combining the X-ray and S- 
Z results and including uncertainties due to the S-Z 
and X-ray normalizations and central temperature, 
we find i/o(go = 5) = 781^8 kms"^Mpc-^ 

There are several additional contributions to the 



uncertainties in yp and Hq. It is possible that the 
cluster gas distribution is aspherical. We use the X- 
ray morphologies of a sample of clusters to estimate 
this contribution to the uncertainty in Hq determined 
from a single cluster. If the cluster exhibits a signifi- 
cant peculiar velocity there is an additional kinematic 
component to the S-Z effect. We use theoretical es- 
timates of cluster peculiar velocities to constrain the 
contribution of the kinematic effect to the uncertainty 
in 2/0 and Hq. An IC gas that contains clumps could 
bias the value of H^. The ROSAT/PSPC + GINGA 
spectrum is used to constrain the contribution that 
clumped IC gas could make to the value of Hq. Fi- 
nally, it is possible that the S-Z effect is astrophysi- 
cally confused. We use theoretical and observational 
results to limit the uncertainty associated with astro- 
physical confusion, including primordial anisotropies 
of the CMB. The contributions to the uncertainty in 
yo and Hf) due to each of these sources are listed in 
Tables ^ and ^. Including these additional sources of 
uncertainty in the analysis which adopts the ASCA 
thermal structure, we determine yo = 3.07J^'gQ x 10^* 
and iJo(go = 5) = 78^40 kms-^Mpc^^ 

There is an important distinction to be made be- 
tween the uncertainties listed in Table ^ those that 
are essentially statistical and those that have a signif- 
icant systematic component. Because we have no in- 
formation about the extent of the IC gas along the line 
of sight, the IC gas is assumed to be spherically sym- 
metric. We must average over a sample of clusters to 
reduce the contribution to the uncertainty due to as- 
phcricity of the IC gas. The uncertainties due to clus- 
ter peculiar velocities and primary CMB anisotropies 
are purely statistical and smaller than that due to 
cluster asphericity. These uncertainties, like the sta- 
tistical uncertainties in the measured data, will be 
reduced for Hq determined from a sample of clusters. 

A more serious obstacle to the determination of Hq 
is the presense of systematic errors which can con- 
tribute to the results for each of the members of a 
sample of clusters. If the IC gas has the thermal struc- 
ture indicated by the ASCA results for A2163, assum- 
ing the gas to be isothermal results in the underesti- 
mation of Hq by ~ 30%. To eliminate this source of 
uncertainty, the thermal structure of each cluster in 
the sample will have to be determined. All IC gas 
will be clumped to some degree; this biases the re- 
sults for Hq to higher values. The presence of clumps 
could be better constrained by high resolution obser- 
vations of the X-ray surface brightness, and improved 
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X-ray spectra. Astropysical confusion of the S-Z ef- 
fect by compact cluster radio sources can be virtu- 
ally eliminated by observations with mm-wavelength 
interferometers. Uncertainty in the mm-wavelength 
flux of the planetary calibrators contibutes ^ ±10% 
uncertainty to the determination of Hq. The CO- 
BRAS/SAMBA satellite mission may be capable of 
relating planetary (and some S-Z) fluxes directly to 
the dipole anisotropy, which is known to a precision 
of ±.2% ( [Fixsen et al. 1996| ). 

To determine an accurate value of Hq using the 
S-Z effect, we must reduce the potential sources of 
systematic error to below the level of the statistical 
uncertainties when averaged over a reasonable sample 
of clusters. Spatially resolved X-ray spectra obtained 
by the next generation of X-ray satellites will be es- 
sential to achieve this goal. Using the S-Z effect, it 
should then be possible to determine Hq to an accu- 
racy of < 10% with a sample of ~ 25 clusters. 
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A. Appendix: Residual Common-mode Sig- 
nal 

In this appendix, we discuss the removal of the 
residual common-mode signal which is not removed in 
hardware due to the mismatch in the responsivities of 
the detectors in a difference. A mismatch in respon- 
sivities of ~ 1% can contribute significant low fre- 
quency noise to a detector difference. Changes in the 
background loading of the detectors over the course 
of an evening can lead to responsivity mismatches of 
this size. Removing the residual common-mode sig- 
nal typically reduces the noise determined from the 
distribution of scan amplitudes by ~ 10 — 20% with 
no significant effect on the fit amplitudes. 

For each difference, we determine the average 
common-mode correlation coefficients (ak) by aver- 



aging akj over the scans in a given period of observa- 
tion; they are typically of order unity. Because of the 
difference in electronic gains between the single and 
differential signals, the amplitude of the source in the 
average single detector signal is smaller than in the 
differential signal by a factor Gg/Gd — 1/120. If Sji 
and dkji were completely correlated we would remove 
a{Gs/Gd) ~ 1% of the differential signal. 

In a drift scan across a stationary source, the dif- 
ference of the signals from two perfectly matched de- 
tectors is orthogonal to their sum. The same is true 
for the signals (dai, ii23, d64j and <456) when com- 
pared to the average of all six single detector sig- 
nals. The fraction of average single detector sig- 
nal to which the differences are sensitive is propor- 
tional to the fractional mismatch of the detector dif- 
ferences, which is ~ a{Gs/Gd)- We estimate the 
amount of signal removed from the differential sig- 
nals as the product of the size of the subtracted sig- 
nal and its coupling to the detector differences to be 
Ay/y « [a{Gs/Gd)f - IQ-*. Subtraction of the 
residual common-mode signal has a negligible effect 
on the determination of the source amplitudes. 
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